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The competition between quantum and classical magnetization plateaus of 5 = 1/2 frus- 
trated Heisenberg chains with modified exchange couplings is investigated. The conventional 
5 = 1/2 frustrated Heisenberg chain is known to exhibit a 3-fold degenerate ][] -type classical 
plateau at 1/3 of the saturation magnetization accompanied by the spontaneous transla- 
tional symmetry breakdown. The stability of this plateau phase against period 3 exchange 
modulation which favors the •-• ^-type quantum plateau state (•-• = singlet dimer) is studied 
by bosonization, renormalization group and numerical diagonalization methods. The ground 
state phase diagram and the spin configuration in each phase are numerically determined. The 
translationally invariant Valence Bond Solid-type model with 4-spin and third neighbor inter- 
actions, which has the exact •-• |-type quantum plateau state, is also presented. The phase 
transition to the classical ||t -type ground state is also observed by varying the strength of 
4-spin and third neighbor interactions. The relation between these two types of models with 
quantum plateau states is discussed. 

KEYWORDS: frustrated Heisenberg chain, magnetization plateau, period 3 modulation, 4-spin interaction, 
Z2 symmetry breakdown, bosonization, numerical diagonalization 
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1. Introduction 

The phenomenon of magnetization plateau has been 
extensively studied as a macroscopic manifestation of the 
essentially quantum effect in which the magnetization M 
is quantized at fractional values of the saturation mag- 
netization Ms in low dimensional magnetism. 

There are at least two types of plateau states that can 
occur for magnetization m = M/M^ = 1/3 in 5' = 1/2 
chain models. One is the "classical state" of Fig. 1(a), 
which occurs as an adiabatic continuation from the Ising 
limit. The other is the •-• | •-• | ... spin configuration 
shown in Fig. 1(b) where •-• stands for the singlet pair. 
This plateau is essentially of quantum origin, driven by 
the tendency of neighboring antiferromagnetically cou- 
pled spins to form entangled singlet states. Since these 
are, in general, only simple caricatures of the ground 
states, we might ask whether they really correspond to 
distinct phases. The ground state symmetries are essen- 
tially the same in both phases. The unit cell consists of 3 
sites, as required by the Oshikawa-Yamanaka- Affleck the- 
orem,^ and there are reflection symmetries about every 
third site and every third link in both cases. To see that, 
nonetheless, that they really are distinct phases note that 
in the classical state the local spin expectation value (Sf) 
is larger on 2 out of 3 sites in a unit cell. On the other 
hand, in the quantum state, it is smaller on 2 out of 3 
sites, the sites where the dimers are drawn. Clearly this 
feature implies distinct phases. 

Recently, Okunishi and Tonegawa*"'^ and Tonegawa 
and coworkers^ investigated the magnetization process 
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Fig. 1. Two possible spin configurations in the m = 1/3 plateau 
state. The classical configuration, (a), is actually observed in the 
model (1) and the quantum configuration, (b), in the model (3) 
for a > 0, with period 3 exchange modulation and also in the 
translationally invariant model (2). 



of the strongly frustrated Heisenberg chain, 

L 

n - JJ2[S^S^+l+5S,S,+2], (1) 

where Si is the spin-1/2 operator on the i-th site. In 
what follows, we take J = 1 to fix the energy unit and as- 
sume the periodic boundary condition Sl+i = Si unless 
specifically mentioned. These authors pointed out that 
this model has a classical plateau at m = M/Mg = 1/3 
for (5 > 5c = 0.487.9 

On the other hand, if we allow for the 4-spin interac- 
tion, we can also construct a model which has exactly the 
quantum plateau by the Valence-Bond-Sohd"^' ^^-type 
construction as follows. 



n 



1 

24 



Y.Tl(Tl-2) 



(2) 



1=1 



where Ti = Si + Si+i + Si+2 + ^'^+3. Each term of this 
Hamiltonian is non-negative definite and projects out the 
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4-spin singlet (T = 0) and triplet states (T = 1) where 
T(T+ 1) = T - . Therefore if we arrange the dimers to sit 
on every third bond as shown in Fig. 1(b), we obtain a 
ground state of this Hamiltonian, since the sum of any 4 
neighboring spins has no projection onto T = 2 states in 
this configuration. 

In the plateau state, one third of the spins remain 
undimerized and they can be in an arbitrary state in 
zero field. These spins are essentially paramagnetic and 
get polarized by an infinitesimal field to yield the 1/3 
plateau state. Although this ground state, at zero field, 
is macroscopically 2^/^-fold degenerate, this degeneracy 
is lifted by any small perturbation. For example, an ad- 
ditional small antiferromagnetic third neighbor antifer- 
romagnetic interaction leads to the Tomonaga-Luttinger 
liquid ground state in the absence of a magnetic field. 
These weakly coupled spin are easily magnetized up to 
the m = 1/3 plateau state. Conversely, a small ferromag- 
netic third neighbor interaction leads to ferrimagnetism 
(i.e. a moment of m = 1/3 at zero field). 

It is also interesting to interpret these two types of 
plateau states in terms of the Jordan- Wigner fermions. 
In the fermionic language, the classical plateau state cor- 
responds to the usual site-centered CDW phase in which 
particles sit on two successive sites and every third site 
is empty, while the quantum plateau phase is a sort of a 
bond-centered CDW or 'bond order wave' (BOW) state 
whc;re the fermion is not localized on a single site but 
tunnels back and forth between a pair of sites and the 
third site is filled. In this context, the distinction between 
classical and quantum plateau states is reminiscent of 
the distinction between CDW and BOW states in iten- 
erant electronic system both of which seem to occur in 
the phase diagram of the half-filled extended Hubbard 
model in one dimension and related models. ^° Consid- 
ering the difference of the nature of these two types of 
plateau states, it is an interesting issue how these two 
plateau states transform into each other. 

The quantum plateau state was actually first 
discovered, over a decade ago in the ferromagnetic- 
ferromagnetic-antiferromagetic (F-F-AF) chain 
model.^'^'^'^ This is a special case of a more gen- 
eral class of models with period 3 space inversion 
invariant exchange modulation: 

L/3 

n = 5^ [(1 - a) (Ss^-i 

1=1 

+ {l + a)S3i+iS3i+2]. (3) 

For a > 1 we obtain the F-F-AF model. The dimers now 
form uniquely on the antiferromagetic links. In this case, 
of course, there is no spontaneous breaking of transla- 
tional symmetry since it is already explicitly broken by 
the Hamiltonian. A more general model with arbitrary 
periodicity is investigated by Cabra and Grynberg.^ 

The competition between quantum and classical 
plateaus is naively realized in the frustrated Heisenberg 
chain with period 3 space inversion invariant exchange 
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modulation described by the Hamiltonian, 

L/3 

W = [(1 - a) {S31-1S31 + S31S31+1) 

1=1 

L 

+ {l + a)S3i+iS3i+2]+6j2SiSi+2. (4) 

i=l 

It is obvious that this model shows a quantum plateau 
for large enough positive a and a classical plateau for 
a = and large enough S. 

This paper is organized as follows. In the next section, 
the low energy effective boson theory for the m = 1/3 
plateau state of the translationally invariant frustrated 
Heisenberg model is presented. The effect of the period 3 
exchange modulation is analyzed using the renormaliza- 
tion group method in §3. The quantitative ground state 
phase diagram at m = 1/3 is obtained by the numerical 
diagonalization method in §4. The nature of each phase is 
also discussed. In §5, the VBS-type model with 4-spin in- 
teraction which has a m = 1/3 plateau is discussed. The 
phase transition to the classical plateau state is investi- 
gated varying the strength of the 4-spin and the third 
neighbor exchange interaction. The physical interpreta- 
tion of the mechanism of the quantum plateau and the 
relationship to the frustrated period 3 model is given on 
the basis of the decoupling approximation. The possible 
universality classes of the transition between quantum 
and classical plateaus in translationally invariant mod- 
els are analysed. In the last section, we summarize our 
results. 

2. Bosonization for Translationally Invariant 

Models 

In order to observe the low energy properties of 
the above models semi-quantitatively, we employ the 
bosonization method^'^' which is generally powerful for 
the description of the one-dimensional spin chains. As 
usual, we first apply the exact Jordan- Wigner transfor- 
mation to lattice fermions, tjjj, then write the low energy 
degrees of freedom in terms of left and right-moving con- 
tinuum fermions with a Fermi wave-vector, kpa = 7r/3. 
(a is the lattice spacing). 

V^, « e'^^/^niaj) + e-'^'^^'^Maj)- (5) 

We see that under the symmetry transformation of trans- 
lation by one site: 

^R^e'^/^^R, iPL^e-"/^iPL. (6) 

This symmetry forbids the usual umklapp term, cx 
(V'1{)^(V'l)^) or more accurately and in accord with Fermi 
statistics: a VH^xV'nV'i^xV'i'- However the "triple umk- 
lapp" term oc is allowed by all symmetries 
and is expected to occur in the low energy effective 
Hamiltonian for the uniform frustrated Heisenberg chain 
(1) and (2) at m = 1/3. Upon bosonizing, we obtain the 
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sine- Gordon model, 



dx 



2/3 



COS 



1 

K J a 

(7) 

where is a boson field, defined on a circle, </> <-> (f)+\/2-K, 
n is the momentum density field conjugate to 0, and K, 
is the Tomonaga-Luttinger liquid parameter. The spin 
wave velocity is set equal to unity. A similar expression 
is also derived in ref. 15. The last term in the Hamil- 
tonian, arising from triple umklapp processes, is rele- 
vant for K < Kc = 4/9. If this term is relevant, the 
Zs-symmetry is spontaneously broken and the ground 
state is 3- fold degenerate as </> = 0)±^^ for 2/3 > 
and = -'-372 2/3 < as shown in Fig. 2. On the 
other hand, if the cos 3-\/2(?!)-term is irrelevant, the ground 
state is the gapless Tomonaga-Luttinger liquid with no 
degeneracy. 

The spin operator Sj on the j-th site and the dimer 
operator on the j-th bond are represented by the boson 
field as, 



4,5 



6 



J+1 



-Ccos (2irj/3 + V2(l)'j 

'Mj - 1) 



const. -|- A cos 



V2^ 



(8) 



where C and A are positive constants of the order of 
unity. Therefore these three ground states correspond to 
the three degenerate m -type and •-• "f-type configura- 
tions for 2/3 > and 2/3 < as shown in Fig. 1(a) and (b), 
respectively. Note that for > 0, the (l) = ground state 
has two larger values of (S'|) (for j = 31 ± 1) and one 
smaller value of (5*1) and also two more negative values 
of {Sj ■ Sj+i) (for j = 31 and 31 — 1) and one less negative 
value. This is precisely what is expected for the classical 
plateau state of Fig. 1(a): {Sj ■ Sj+i) is more negative 
on the bonds joining spins ordered in opposite directions. 
Conversely, for ys < 0, the </> = tt/v^ ground state has 
only one larger value of (S'j) (for j = 31) and one more 
negative value of {Sj ■ Sj+i) (for j = 3/ + 1), correspond- 
ing to the quantum plateau state of Fig. 1 (b) . Although 
the sign of 2/3 cannot be analytically determined, we may 
conclude 2/3 > for the Hamiltonian (1) for large enough 
S, because ttitti ••• structure is actually observed by 
Okunishi and Tonegawa.^ On the other hand, the case 



1-i-a 1-a 1-a 1-i-a 1-a 



(b)l 
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l+a 1-a 1-a 1-i-a 1-a 
or 

t t + t t 

1-Hx 1-a 1-a 1-i-a 1-a 



t t I t t I 



1-i-a 1-a 1-a 1-i-a 1-a 

Fig. 4. The spin configuration in each phase: (a) yi > j/ic, (b) < 
yi < yic, {c)yi < 0. The pairs of spins connected by thick lines 
form singlet pairs. The three spins in rectangles form doublets. 



2/3 < corresponds to the model (2) which we argued 
exhibits the quantum plateau. 

3. Period-3 Exchange Modulation 

3.1 Bosonized Hamiltonian 

Now we consider the effect of period 3 exchange mod- 
ulation in the models (3) and (4). Since these models 
exhibit the classical plateau at a = and large enough 
6, we assume that the parameter, 1/3 > 0. Exchange mod- 
ulation is included in the bosonized Hamiltonian by the 
perturbation,^' ^' ^ 



Hi = J^^ , / dx cos 
27ra^ J 



(9) 



where yi = AnaaA/S oc a in terms of the boson field. 
This term reduces the 3-fold degeneracy to 2-fold degen- 
eracy for 2/1 > as long as \yi \ is small as depicted in Fig. 
3. For 2/1 < 0, this term totally removes the degeneracy. 
In the former case 2/1 > 0(a > 0), the ground states are 

= ± ' 



=atan — 



Qya-yi 



which are close to ' 



3V2 



VS"""" V y ^ya+yi 
for small \yi\. Here the branch of arctangent is taken as 
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atana; € (0,7r]. The spin configurations of tlicsc two de- 
generate states for small |yi| are nearly ittitt ••■ and 
TiTTiT ••• according to (8) as shown in Fig. 4(b). This 
implies that the Zi translational symmetry is sponta- 
neously broken in this state. If y\ exceeds yic = 9y3, 
these two ground states merge to a nondegenerate ground 
state cf) = -K / \/2- This state is the quantum plateau state 
with •-• I •-• 1" ... spin configuration shown in Fig. 4(a) 
according to the expression (8) which is essentially the 
same as the ground state of (3) with a > 

In the latter case yi < 0(a < 0), the ground state 
is = 0. According to the expression (8), this phase is 
the classical plateau state with spin configuration ttitti 
... as shown in Fig. 4(c). In this phase, the three spins 
coupled via stronger exchange interaction (1 — a) form 
antiferromagnetic trimers whose ground state is the S = 
1/2 doublet tit and |ti- bi the plateau state, the whole 
system is covered by the tit trimers resulting in the 
classical plateau state with the symmetries of the simple 
TTitn... state. 

The transition between the phase with spontaneously 
broken Z2 symmcrtry and the; •-• | •-• | ... -phase is ex- 
pected to be of the Ising type from symmetry consider- 
ations. On the other hand, the transition at yi = is 
a first order transition if 7/3 is relevant [5 > Sc) because 
the value of jumps at the transition point. If ys is ir- 
relevant {S < Sc), the transition at yi = is a Gaussian 
transition because the relevant nonlinear term is absent 
in this case. 



3.2 Renormalization Group Treatment 

For yi < 0, the cos •\/2(/)-term and cos3-\/20-term do 
not compete. Therefore the ground state is obviously 
unique (f> = and gapful as long as cos ■\/2</'-terni is rel- 
evant. However, for yi > 0, these two terms compete with 
each other and drive the Ising type phase transition. To 
describe this competition, it is necessary to resort to the 
renormalization group analysis. Up to the second order 
in yi and ys, we find the renormalization group equation, 

d 1 _ 9 2 



dlK 
dys 

dl 
dyi 

dl 



= (2-^)2/3, 
= (2-f)2/i, 



(10) 

(11) 
(12) 



taking into account the scaling dimensions and the coef- 
ficients of the operator product expansion(OPE) of the 
operators {dx4>)'^, cos-\/2(/> and cos3\/2(/). Confining our- 
selves to the region of small yi cx a, we have neglected 
the terms of the order of y\ . Because we are interested in 
the behavior near the plateau-non-plateau critical point, 
we set K = |(1 -|- ^yo) and keep the lowest order terms 
in yo- 

dyo 
dl 

dm 
dl 

dyi 



= -2/3, 



-yoys, 



dl ^99 



(13) 
(14) 
(15) 



Similar set of equations are obtained and analyzed by 
Kadanoff^^ and Kitazawa and Nomura^^ for the double 
sine-Gordon model with cos a/20 and cos 2y/2(j)-tenns. In 
the present case, the coupling between and yi does not 
appear in the rhs of (15), because the operator cos\/2(/) 
does not appear in the operator product expansion of 
cos\/20 and cos3\/20- The first two equations are the 
well-known Brezinskii-Kosterlitz-Thouless renormaliza- 
tion group equations. The solution which describes the 
plateau state is given by 



2 2 

ys-vo 



t > 0, (16) 
yo = -\/itan[\/i(Z-/o)], (17) 



= r/T,, ,,1 - (18) 

cos[vHt — tojj 

where t measures the distance from the critical separatrix 
Vz~yo = 0. In the present case t oc S—Sc- The bare values 
of yi, yo and y3 are denoted by yio ^ a and y3o,yoo ^ 
Oil). 

The integration constant /g is determined by the initial 
condition. 



2/30 



cos[\/tZo] 



0{l) 



(19) 



This implies cos\/t?o ^ 0{Vi) and Iq ~ 7r/2\/t + 0(1) 
for t « 1. Equation (15) can be integrated using these 
solutions to give 

ln(yi/yio) + l Hys/yso)- (20) 

for t << 1. The phase transition takes place by the com- 
petition between yi and y^. Therefore the phase bound- 
ary ac{S) can be estimated by setting yi ^ y^ ^ 0{1) 
after enough steps of renormalization, that is at some 
value of I. Since yi{l) grows exponentially with I, from 
Eq. (20) but ysil) grows more slowly, from Eq. (18), be- 
coming 0(1) at / « ir/Vt, we see that yi(7r/-\/i) must be 
0(1) on the phase boundary. This implies: 



(21) 



or 



Qc « C'e-^/^^, (22) 
where C and C are constants of the order of unity. 

3.3 Numerical Diagonalization Studies 

3.3.1 Ising Transition (S > Sc,a = a^S) > 0) 

To confirm the phase diagram quantitatively, we carry 
out the numerical diagonalization studies of finite size 
chains. 

As discussed earlier, the transition between the Z2- 

symmetry broken intermediate phase and nondegener- 
ate phase for large a belongs to the Ising universality 
class. The energy gap between the ground state and 
the lowest excited state with M = Ms/ 3 should be- 
have as A (X const, in the non-degenerate phase and 
oc exp(— const, x L) in the Z2-broken phase. (The expo- 
nentially small energy gap for a finite system results from 
tunnelling between the two degenerate ground states.) 
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On the other hand, at the critical point, A oc 1/L. 
This impHes that the phase boundary can be success- 
fully determined by the phenomenological renormaliza- 
tion group (PRG) method. The finite size critical point 
ac{5\ L,L + 6) can be determined by, 

LA{ac{6;L,L + 6),L) 

= {L + 6)A{ac{5;L,L + 6),L + 6) (23) 

for fixed S. We take the system size L as multiples of 
6 {L = 12,18,24 and 30). The critical point adS) in 
the thermodynamic limit is determined by extrapolating 
ac{S;L,L + 6) to L — > oo as shown in Fig. 5 assuming 
the formula,i«'i9 



ad5;L,L + 6) 



const 
(L + 3)3' 



(24) 



as plotted in Fig. 6. Here we have taken into account 
that the present transition is expected to belong to Ising 
universality class from symmetry consideration. 

The exponent v of the energy gap A defined by A ~ 
(ac — a)" is also estimated by the PRG method as, 

ln((L + 6)/L) 



u{L,L + 6) = 



In 



{L + 6)A' {ac{5; L, L + 6), L + 6) 
LA'{ad6;L,L + 6),L) 



(25) 

Here A' is the derivative of the energy gap with respect 
to a at fixed L. Extrapolation to L — > oo is carried out 
assuming the formula,^®' 

const 

u{L,L + 6) ~ (26) 

as plotted in Fig. 7. Although the size dependence is 
rather strong, we have verified that the value of f is close 
to the Ising value {v = 1) as long as the value of 5 is away 
from Sc- 

In Fig. 8, \nac{S) is plotted against l/VfJ — with 
5c = 0.487. This plot shows that the 5-dependence of 
the; numerically obtained value of ac is well described 
by the formula exp(— const./-y/|5 — Sc\) derived by the 
renormalization group method in the preceding section. 

3.3.2 First Order Transition (6 > Sc, a ^ 0) 

For S > Sc, the numerically obtained ground state 
energy of the finite size system shows an almost level 
crossing-like behavior at a = as shown in Fig. 9 sug- 
gesting the first order transition. Although this is not 
an exact level crossing for finite L, the energy eigen- 
values per site of two low lying states are almost size- 
independent near the transition point. Therefore we can 
interpret these two states as two macroscopically distinct 
states which are adiabatically connected to the ground 
states for a > and a < and the first order transition 
takes place between these two states. This is consistent 
with the bosonization argument in the preceding section. 

3.3.3 Gaussian Transition (S < ScCt = 0) 

For S < Sc, no level crossing between two macroscopi- 
cally distinct states takes place suggesting the continuous 
transition as shown in Fig. 10. We have studied the tt- 
twist boundary condition which corresponds to changing 




0.0002 l/(L+3)^ 0.0004 



Fig. 5. Extrapolation procedure of the finite size critical points 
ac(5:L,L + 6) to L ^ oo. The dominant size dependence is 
assumed to be -|- 3)^. The values of 5 are 1, 0.95, 0.9, 0.85, 
0.8, 0.75, 0.7, 0.65, 0.6 and 0.55 from top to bottom. 



a 
0.8 



0.4 



non-plateaus^ 



5=0.487 




Fig. 6. The ground state phase diagram at m = 1/3. 




0.002 



0-004i/(L+3)' 



Fig. 7. Extrapolation procedure of the finite size critical exponent 
v{L, L + Q) to L — » oo. The dominant size dependence is assumed 
to be 1/(L 4-3)2. 



the sign of the [S+S^ + 3+ Si), {St_^Sl + 5'+S'^_J 
and {8^82 + 82 S^) terms in the Hamtilonian, find- 
ing that the space inversion parity of the ground state 
reverses at a = 0. This is a clear indication of the Gaus- 
sian transition driven by the term (9) as discussed by 
Kitazawa^° and Kitazawa and Nomura. We have also 
numerically checked that the conformal charge c is close 
to unity along the line a = and S < Sc &s shown in 
Fig. 11. Therefore this transition is confirmed to be the 
Gaussian transition. 
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Fig. 8. The critical value of In a plotted against l/\/S — Sc 



-0.28 
E/L 



-0.285 



-0.2S 



5=0.7 



.01 



0.01 



Fig. 9. The lowest two energy eigenvalues per site at 5 = 0.7 for 
L = 30(open circles) L = 24(open squares) and L = 18 (filled 
squares). The data points for L = 24 and L = 18 are almost 
covered by those for L = 30. 
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5=0.3 a ■ I I ■ ■ ■ 

ggeQQOOOOOO 
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000000000399 



-0.01 



0.01 



Fig. 11. The ^-dependence of the numerically estimated confor- 
mal charge for a = 0. Finite values of c for 5 > 5c are numerical 
artifacts. 
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10 
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Fig. 12. The spatial variation of (5?) in the ground state of the 
finite chain with L = 30 for {5, a) = (0.8, — 0.1)(open circles), 
(<5, a) = (0.8, 0.1) (filled squares) and (0.8, 0.1) (open squares). 



Local field H 



loc 



_1-a _1+a_1-a_1-a _ 
• • • • • 

= L 1 2 3 4 



Fig. 13. Local field applied to the end spin. 



Fig. 10. The lowest two energy eigenvalues per site at <5 = 0.3 for 
L = 30(open circles) L = 24(open squares) and L = 18 (filled 
squares). 

3.3.4 Nature of the Z2-broken phase 

To verify the nature of each phase described in §2, 
wc numcricany calculate the spin configuration in the 
plateau state. 

In Fig. 12, we plot the expectation value {Sf) for 
((5, a) = (0.8,0.4), (0.8,0.1) and (0.8,-0.1) against i for 
L = 30. The spin configuration for {5, a) = (0.8,-0.1) 
shows a clear TTiTTi ••• structure. On the other hand, 
those for {8, a) = (0.8,0.1) and (0.8,0.4) look quite sim- 
ilar indicating the •-• f f ... structure at first sight. 
However, even in the Z2-symmetry broken phase, the 
ground state of the finite size; system does not break 
the symmetry but is the antisymmetric linear combi- 
nation of the symmetry broken states littiTT ••• ) and 
ITiTTit ••• )• Therefore it is difficult to distinguish these 
two phases by naive numerical calculation for the finite 



size systems. 

In order to single out the ^2-symmetry broken phase, 
we apply the local upward field magnetic field Hioc on 
the i = 1 site as shown in Fig. 13. The calculated ex- 
pectation values of (S'f ) are shown in Fig. 14. The Z2- 
symmetry broken state ITiTTit ••• ) is clearly singled out 
for {6, a) = (0.8,0.1) with weak local field Fioc = 0.02. 
On the other hand, the •-• f •-• f ... spin configuration 
for {6, a) = (0.8,0.4) and the TTiTTi ■•• spin configura- 
tion for (^, a) = (0.8,-0.1) are stable against local field 
except for the tiny polarization. It should be noted that 
no spatial periodicity is imposed by the local magnetic 
field in the present calculation. 

In Figs. 15 and 16, the local field dependence of the 
magnetization on the central sites i = L/2 + 1 and 
L/2 + 2, which are furthest from the position of the lo- 
cal field, are presented for {S,a) =(0.8,0.1) and (0.8,0.4), 
respectively. For {6, a) = (0.8,0.1), the magnetization 
saturates for small local field and this saturation value is 
insensitive to the system size. This implies these are the 
bulk values of local magnetization. On the other hand. 
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Fig. 14. The spatial variation of (Sf) in the ground state of the 
finite chain with L = 30 for {S,a) = (0.8, — 0.1)(open circles), 
(5, a) = (0.8, 0.1)(filled squares) and (0.8, 0.4) (open squares) 
with small upward local field ifioc = 0-02 at i = 1 site. 
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Fig. 15. The local field dependence of (Sf ) on the central sites 
for (S, a) = (0.8, 0.1). The filled symbols are ^^'^ open 

symbols are {^^ 1^/2+2)' 
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Fig. 16. The local field dependence of {Sf) on the central sites 
for (5, a) = (0.8, 0.4). The filled symbols are (^\i2-\.\^ ^'^'^ open 
symbols are (j^ 1^12-12) ■ 



for {5, a) = (0.8,0.4), the magnetization does not satu- 
rate for small local field. Fig. 17 shows the system size 

dependence of the initial slope of the iJioc-dependence 



of is 



'L/2+1 



and ( >5'£^2+2 / ^'^^ 



(0.8,0.4). It is 

evident that the slope decreases exponentially with the 
system size. This confirms that the bulk ground state is 
actually the |»-» | •-• | ... ) state in this phase. 




Fig. 17. The size dependence of the initial slope 
d{S?) /dHioc\Hioo=o for (5, a) = (0.8,0.4). The symbols 
for i = L/2 + 1 and i = L/2 + 2 overlap within the size of the 
symbols. 



4. VBS-type Model 

Another way to have a quantum plateau is to reverse 

the sign of the cos term in the bosonized language. 

In this case, the quantum plateau is induced by the spon- 
taneous symmetry breaking in the uniform Hamiltonian. 
However, within the two-spin exchange interaction we 
do not find the appropriate model in which quantum 
plateau is realized. Instead, as explained in section 2, 
we can construct a translationally invariant VBS-type 
model (2) with 4-spin interaction which exactly shows 
the 1/3 quantum plateau with spontaneously broken 
translational invariance. 

In terms of the original spin operators Si, the Hamil- 
tonian (2) is rewritten as. 



n = 

+ 
+ 
+ 



1 ^ 

- ^[30Si5i+i + 20SiSi+2 + lOSiSi+3] 



24 

8[(S,Si+i)(Si+2iS'i+3) 

{SiSi+2){Si+lSi+3) 
{SiSi+3){Si+iSi+2)] 



(27) 



except for the constant term. 

In order to make clear that this quantum plateau state 
is a distinct phase from the classical plateau found in 
model (1), we consider the Hamiltonian which continu- 
ously interpolate between the Hamiltonians (1) and (27) 
as follows. 



n = 



1 

— J2[30SiSi+i + 20SiSi+2 + lOeSiSi+a] 



i=l 



89[{SiSi+i){Si+2Si+3) 

("S^i 5^1+2 ) ('S'i+ 1 ) 
('S'iS^i+3) (S'i+i S'i-|_2)] 



(28) 



and investigate the ground state transition with respect 
to 6. In the present case, however, the phases on both 
sides of the transition are phases with spontaneously bro- 
ken Z3 symmetry. Therefore the lowest energy gap de- 
pends on the system size as A ~ exp(— const, x L). This 
implies that the naive application of the phenomenolog- 
ical renormalization group analysis is not possible. How- 
ever, consider the effect of open boundary conditions, 
for a chain of length L = 6n, in the two phases. In the 
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classical plateau phase we expect that open boundary 
conditions will favour a unique plateau state, raising the 
energy of the other two by an L-independent amount. 
The favoured state is the one drawn in Fig. 18(a) in the 
case n = 1, L = 6. Note that only 2n — 1 energetically 
costly parallel neighboring spins occur in this state. The 
other two classical plateau states have 2n parallel spins. 
On the other hand, consider the quantum plateau states 
with open boundary conditions. We expect that the en- 
ergy will be lowest if no singlet bonds are broken. In this 
case a singlet bond will occur at one end of the system 
and a polarized spin at the other, as indicated in Fig. 
18(b) for L = 6. This allows two degenerate states, re- 
lated to each other by space inversion. For a finite open 
system we expect tunnelling process to mix the two low 
energy quantum plateau states, of Fig. 18(b) leading to 
an exponentially low energy excited state. On the other 
hand, in the classical plateau state the gap should be 
0(1). 

Therefore the energy gap within the m = 1/3-sector 
behaves as A ~ const, in the former phase and A ^ 
exp(— const. i) in the latter phase and PRG analysis is 
possible. The extrapolation procedure of the finite size 
critical point 9c{L,L + 6) determined from the cross- 
ing point of LA for the sizes L and L + 6 is shown in 
Fig. 19. The transition point in the thermodynamic limit 
is 6c — 0.57. Although the precise value of the critical 
point is not accurate enough due to the relatively large 
size dependence of finite size critical point, this analy- 
sis clearly shows that the •-• t and m states are two 
distinct phases even in the uniform chain. 

This transition can be interpreted in terms of the 
bosonized language. At 9 = Oc, the coefficient j/3 van- 
ishes. Therefore, the higher order umklapp term 



H 



ye 

27ra2 



da; cos 



(29) 



comes into play. This term is irrelevant or relevant de- 
pending on whether i('>l/9orii'<l/9. If this term 

is irrelevant, the transition is expected to be the second 
order Gaussian transition with conformal charge c = 1. 
If this term is relevant and ye < 0, the ground state is 6- 
fold degenerate at the transition point. This corresponds 
to coexisting classical and quantum plateaus. Turning 
on j/3, with either sign, either the classical or quantum 
plateau states arc favored. In this case, the transition is 
a first order transition and level crossing behavior is ex- 
pected. If this term is relevant and ye > 0, the 2/3 and 
2/6 terms compete near the point ys = and two Ising 
type transitions should take place as discussed in the 
competition between cos and cos 2\/2(?!) terms in the 
Ashkin-Teller model.i^-i'^ 

Numerically we find only a single transition and the 
smallest energy gap A with periodic boundary conditions 
scales as LA ^ const, near the critical point. Wo also 
find the conformal charge c ~ 1.0 extrapolating from 
L = 12, 18 and 24. Therefore we conclude that the first 
one of the above three senarios is valid in the present case 
and the transition is a Gaussian transition. In order to 
confirm that ye is irrelevant and j/3 is relevant near the 
transition point, we determine K using the relationship 



(a) 



t I t t ; t 



t 



t 



t 



Fig. 18. The ground state spin configurations of the Hamiltonian 
(28) with open boundary conditions and 6 sites. 




2/(U3) 



Fig. 19. Extrapolation procedure of the finite size critical point 
ec{L,L + 6) to L -> 00. 



between the excitation spectrum and K. Here we employ 
the relation which contains the lowest excitation energy 
of the mode with wave number k = 2^/3, 

K _ LA{k = 2Tr/3,M = M,/3) 

y ^ 



(30) 



(31) 



2ttVs 

where Vs is the spin wave velocity given by, 

_ LA{k = 2tt/L, M = Msl3) 
2^ • 
Extrapolating the value of K{L) estimated for finite size 
L to L ^ 00 as K{L) ~ if + C/L^ from L = 12, 18 
and 24, we find K ~ 0.368 which is between 4/9 and 
1/9. Therefore 1/3 is relevant and ye is irrelevant and the 
condition for the Gaussian transition is fulfilled. 

To understand the physical origin of the •-• t phase in 
the 4-spin Hamiltonian, we decouple the 4-spin terms in 
(27) as, 

{S,S,){SkSi) ^ {{S,S,)) (SkSi) 

+ {S,S,){{SkSi))-{{S,S,)){{SkSi)) (32) 

Assuming the complete •-• f •-• f ... state, we take 

3 

{{S 31+1 S 31+2)) = -7 



{(SsiSsi+s)) 
{(SiSj)) 



1 

4 

other terms 



(33) 



Then we have 



L/3 



= Y^iy^iS3l-lS3l + S3lS3l+l) 
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Fig. 20. The spatial variation of (S?) in the ground state of the 
decoupled model with L = 30 with small upward local field 
ifloc = 0.02 at i = 1 site. 
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Fig. 21. The local field dependence of {Sf) on the central sites 
in the ground state of the decoupled model . The filled 

symbols are ^'S'£^2+i) ^"^"^ open symbols are ^'S'£^2+2)' 



+ 32Ssi+iSsi+2 

+ 20(S3;_iS3i+i + S3;S3;+2 + S3I-2S31) 
+ 10{S3l-2S3l+i + S3I-1S31+2) 

+ 4S3;S3;+3] (34) 

except for the constant terms. If the last term is ne- 
glected, this Hamiltonian is equivalent to the Hamilto- 
nian (4) with a = i ~ 0.143 and 5 = | ~ 0.714 which 
belongs to the •-• "f phase in the phase diagram Fig. 6. 
To demonstrate that the last term is not harmful to the 
•-• I •-• I ... structure of the plateau, we carried out the 
numerical diagonalization of the Hamiltonian (34) and 
calculated the local spin configuration in the ground state 
with local field at z = 1. The result is shown in Fig. 20 for 
L = 30. The expectation values of the dimer spins almost 
vanishes and the polarized spins are almost equal to 1/2. 
We also calculated the spin-spin correlation functions 
which we have assumed in Eq. (33) as, 

((S3i+iS3/+2)) = -0.74294967 

((S3(S3;+3)) = 0.24014714 

I {{SiSj)) I < 0.04020856 other terms (35) 

for _L = 30. The results agree with the assumed set of 
values (33) quite well. Therefore we may conclude our 
assumptions (33) are self-consistent and we can describe 
the physical picture of the quantum plateau state in the 



following way. 

The periodic array of dimers induce the effective pe- 
riod 3 exchange modulation via the 4-spin term self- 
consistently, resulting in the quantum magnetization 
plateau with spontaneous Zs-translational symmetry 
breakdown at m = 1/3. This implies the 4-spin inter- 
action can be generally the source of dimer type order in 
low dimensional magnets. 

5. Summary 

The competition between the m = 1/3 classical and 
quantum magnetization plateau states in the 5=1/2 
frustrated Heisenberg chains with space inversion invari- 
ant period 3 exchange modulation and 4-body interac- 
tion is investigated by the renormalization group and 
numerical diagonalization method. 

The conventional S = 1/2 frustrated Heisenberg 
chains is known to exhibit the 3-fold degenerate tit - 
type classical plateau accompanied by the spontaneous 
Z3 translational symmetry breakdown. Two of them 
(littitt ••• ) and ITiTTiT ••• )) turned out to be robust 
against the period 3 exchange modulation which favors 
the •-• I •-• I ... phase up to a critical value of the mod- 
ulation amplitude (•-• = singlet dimer) resulting in the 
Z2 translational symmetry broken phase. Another ||| - 
type state ITTiTTi-") is stabilized for period 3 mod- 
ulation with opposite sign. The transition between the 
•-• I •-• I ... -phase and Z2-broken phase is the Ising 
transition and that between the TTiTTi ••• -phase and 
Z2 broken phase is the first order transition. The spin 
configuration in each phase is numerically verified. 

A translationally invariant Hamiltonian with the ex- 
act •-• ^-typc quantum plateau ground state is also pre- 
sented based on the VBS-type construction. The sponta- 
neous Z3 translational symmetry breakdown takes place 
in the ground state. The phase transition from the •-• f- 
type phase to the -type phase takes place as the 
strength of the 4-spin term and third neighbor interac- 
tion is varied. This implies that these two phases are 
distinct phases. The universality class of this transition 
is discussed in bosonized language. From numerical anal- 
ysis, this transition turned out to be the Guassian tran- 
sition. A physical picture of the quantum plateau phase 
is presented based on the decoupling approximation. 
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